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Abstract 

A general scheme has been proposed to study the critical behaviour of integrable 
interaction-round-a-face models with fixed boundary conditions. It has been shown 
that the boundary crossing symmetry plays an important role in determining the 
surface free energy. The surface specific heat exponent can thus be obtained without 
explicitly solving the reflection equations for the boundary face weights. For the 
restricted SOS L-state models of Andrews, Baxter and Forrester the surface specific 
heat exponent is found to be a s = 2 — [L + l)/4. 







1 Introduction 



Since Sklyanin's work |J integrable models with open boundary conditions have received 
considerable attention (see, e.g., || [| || ^| and references therein). However, despite 
this flurry of activity, much of the exact surface critical behaviour of the various models, 
both in quantum field theory and in statistical mechanics, remains to be derived. Our 
motivation here is to apply the boundary crossing symmetry investigated by Ghoshal and 
Zamolodchikov in terms of quantum field theory || to study surface critical phenomena 
in integrable lattice models in statistical mechanics. 

The particular models of interest are the interaction-round-a-face (IRF) solid-on-solid 
(SOS) models. First among these is Baxter's SOS model [[?[] which is intimately related 
to the eight- vertex model ||. In this SOS model the heights round a given face are 
unrestricted. Under restriction of the heights, this model becomes the L-state restricted 
SOS model of Andrews, Baxter and Forrester (ABF) ||. There has been increasing 
interest in studying these SOS models with open boundary conditions |5], |It], [IT| . In this 
case the integrability governed by the usual star-triangle equation must be supplemented 
by the reflection equations. In [|1IJ the intertwining relation between the vertex and face 
weights at the boundary has been studied and thus, starting from the reflection equations 
of the eight-vertex model, the IRF analogue of the reflection equations and the general 
solutions of the boundary face weights have been formulated using the boundary face- 
vertex intertwining relation. Reflection equations for IRF models have also been written 
down in @, |, [Tl[ . 



In this paper we show that the boundary crossing symmetry of the boundary face 
weights is enough to give the unitarity relation for the surface free energy on a square 
lattice rotated by 45° with fixed boundary conditions. Thus we can extract some surface 
critical behavior without knowing the explicit details of the boundary face weights. The 
underlying models considered here are the unrestricted SOS model and ABF restricted 
SOS models. However, the discussion is applicable to other integrable IRF models. 

In the next section we describe the models with open boundary conditions. Instead 
of solving the SOS reflection equations we assume that the generic solutions of the re- 
flection equations define the boundary face weights, which satisfy the boundary crossing 
symmetry. In section [5] we show how to construct the unrestricted or restricted SOS 
models with fixed boundary conditions from the models with the open boundary condi- 
tions. This is done by taking the special open boundary face weights and introducing 
alternating inhomogeneities such that the square lattice is rotated by 45°. In section |] 
we present the functional relations for whole fusion hierarchies of the models. Ignoring 
the finite-size corrections to the transfer matrices we end up with a unitarity relation 
determining both the bulk and surface free energies. This relation is simplified by the 
boundary crossing symmetry. Then in section [| we study the surface critical behavior 
of the models with fixed boundary conditions. In particular, the surface critical expo- 
nent corresponding to the specific heat does not rely on the details of the boundary face 
weights. As a by-product, using the known scaling relation, we also obtain the correlation 
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length exponent v — (L + l)/4 for the ABF model with odd L. 



2 Models with open boundary conditions 

The square lattice ABF models are restricted solid-on-solid (RSOS) models with L heights 
built on the classical Dynkin diagram. The corresponding unrestricted SOS model 
has been introduced by Baxter in the study of the eight- vertex model J7j. By using the 
intertwiners the eight-vertex model can be transformed into the SOS model, which is 
defined by the following nonzero face weights 
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where w a = aX + w - The height variable a G Z with w and A as free parameters. The 
function $i (u) is a standard elliptic theta functions of nome p, 
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The above face weights satisfy the star-triangle equation 
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along with inversion/unitarity/crossing unitarity relations 
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where p(u) = i? x (A - u)$i(A + u)/$\(X). 

The A L ABF RSOS models follow from the unrestricted SOS model on setting w = 
and A = with a = 1, 2, • • • , L, where L = 3, 4, • • •. In particular, the Ising model is 
described by the A 3 model. 

The integrable boundary face weights are represented by a triangular face with three 
spins interacting round the face [|5], [TO], 
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They satisfy the reflection equation ||, |T^ 
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where generally we may have more arbitrary parameters than one like £. The SOS 
analogue of boundary crossing symmetry is given by 
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This relation is expected by applying the intertwining relation to the boundary crossing 



symmetry of the eight-vertex model flq, 10 



Following [10 we define the row-to-row transfer matrix T(u) with elements 
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where a = {ao, ai, ■ ■ • , a^v} and b = {bo, &i, - • • , &v}- Here are arbitrary parameters 
which play the role of inhomogeneities. The right boundary face weights K_ are given 
by 
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and the left boundary face weights K + are given by 



/\ . ( ^ a 



u 



K a 



-« + A;^- 



Xw b )'d 1 (w c 



So defined, the transfer matrix forms a commuting family 

[T(«) , T(v) ] = . 
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It follows that the above SOS and RSOS models with open boundary conditions formu- 
lated by the boundary face weights K are integrable systems. 



3 ABF models with fixed boundary conditions 

The square lattice rotated by 45° is a natural geometry to investigate surface critical 
phenomena. The ABF model in this geometry has received particular attention W2 
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Figure 1: The rotated square lattice following from the choice fl3.1| ) for the alternating inho- 
mogeneities in the unrotated square lattice. Each face has spectral parameter 2u while the 
boundary faces have spectral parameter u. 




Figure 2: The rotated square lattice with the boundary spins a,b,c,d fixed: \a — b\ = 1 and 
| c — cZ| = 1. The boundary face weights are dropped after taking the free parameter £ = £q- 



Dj. From the perspective of exactly solved vertex models it is known that the rotated 
geometry can be realised by appropriate choices of the inhomogeneities in the unrotated 
lattice [14|, [|. Here we apply the same idea to the ABF model, however the discussion is 



also valid for the unrestricted SOS model. We show that the ABF model on the rotated 
lattice with fixed boundary conditions follows from the ABF model on the unrotated 
lattice with the corresponding open boundary conditions. 

Let us consider the ABF models on the unrotated square lattice with the row-to-row 
transfer matrix denned in ( |2.9|) . For simplicity we consider only the diagonal i^-matrices, 
or alternatively, the boundary face weights (|2.6|) vanish if b ^ c. The particular choice of 
inhomogeneities 

v k = (~l) k u (3.1) 
leads to the desired orientation because of the property 
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To see this it is most instructive to view the situation graphically as was done for the 
vertex models in ||. In this way we arrive at the rotated lattice in Fig. 1, where the face 
weights have been depicted graphically as 
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The boundary face weights represented by the triangular faces contain the free parameter 
£ which, following || can be chosen such that 



K a 







for a — b + 1 or a — b — 1 . 



Thus the boundary face weights can be dropped under normalisation. In this way we 
arrive at the ABF models on the rotated square lattice with fixed boundary conditions, as 
depicted in Fig. 2. This is the special geometry considered in [O. 1 3 1 with the particular 



fixed boundary conditions corresponding to arrow conservation at the boundaries in 
the vertex formulation. Like their vertex model counterparts the procedure shown here 
provides a possibility to exactly solve an integrable SOS or RSOS model on the rotated 
square lattice with fixed boundary conditions by first solving the corresponding model 
on the unrotated square lattice. Thus the coordinate, algebraic or analytic Bethe ansatz 
could play an important role again. 



4 Fusion hierarchies and functional relations 



The fusion procedure for constructing integrable generalisations of the SOS or ABF mod- 
els has been described in [|1|] (see [[L6j for a related work). The fusion of integrable open 
boundary face weights can be done in a similar manner to form fusion hierarchies of new 
integrable models defined via the fused weights [|Tj], [18], || [l(| . The functional relations of 
the fusion hierarchies of the ABF models with diagonal open boundary conditions (b = c 
in Q2.6|) ) have been given in ||. However, the transfer matrix defined in |J is different to 
the one considered here, so it is a worthwhile exercise to derive the functional relations 



for the fused transfer matrix defined in ( |4.1| ) with general boundary fused face weights 



based on ( pT6| ) . The derivation is in the same spirit as in |5|, [L0| and is to apply the fusion 
procedure to the product of two transfer matrices. 

For a square lattice the fusion can be completed separately in the vertical and horizon- 
tal directions. The fused face weights are expressed by W(n,m) ( M ) ^ ne f us i° n 

levels n and m, respectively, in the vertical and horizontal directions [IB]. The fused 
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boundary face weights are given by | a 

It has been shown that the fused face weights satisfy the star-triangle equation and fused 
boundary face weights satisfy the reflection equation |5|, |TU|. Thus the fused models with 
open boundary conditions are still integrable if the corresponding fused transfer matrices 
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They form fusion hierarchies of commuting families of transfer matrices 







(4.2) 

where the fusion level m labels the different families and fusion level m = n = 1 is the 
unfused transfer matrix. The fused face weights vanish when the fusion level m is greater 
than L — 1 only for the ABF model [IB|. So we have L — 1 families in total for the 
ABF model while for the SOS model we have an infinite number of families. In each 
family the corresponding transfer matrices are related to each other and satisfy a group 
of functional relations. To see this let us define 



T (n) = T {m,n)^ + kX) m = 1? 2> . 

T (n) = if n < or m < 
T (0) = I 
Then the functional relations read 
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with no closure condition (n, m = 1, 2, • • ■) for the SOS model but with 
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as the closure condition (m = 1, 2, • • • , L — 1) for the ABF model. The matrix function 
f™{u) is related to the anti-symmetric fusion and is dependent on both the bulk and 
boundary face weights. 

The functional relations can be proved in a similar manner to || l0| . We find the 
matrix function f™{u) to be given by 

f?(u) = f m (u + nX) (4.5) 

r» = nW-iA) P ( M +jA)r, (4.6) 

P\^ U ) j=Q 

where the boundaries contribute the factors uj~(u) and u + (u), which are diagonal matri- 
ces labelled by the heights along the right (r) and left (/) boundaries, 
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Applying the boundary crossing symmetry to the surface face weights we have 

tfi(2A + 2w) 
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Similarly oj + (u) is given by uj (u) with u —>■ —u, £_ — > and r — *■ I. According to the 
spirit presented in [|It], [HJ the crossing unitarity condition 



T(u)T(u + A) = f\u)p(2u)$\(\)/$\(2\) 
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determines both the bulk and surface free energies, where p{2u)'&\{\) j '"&\(2X) has been 
used for the normalisation of the surface free energy. 

5 Surface critical phenomena 

The surface free energy can be obtained by applying the "inversion relation trick" , which 
is known to give the correct bulk free energy of the eight- vertex model [ 2(J . The bulk and 
surface free energies must both satisfy the unitarity relation ( |4.10| ) and can be separated 
from one another. In the following we consider only diagonal A-matrices. 



Let A (it) = Ab(u)A s (u) be the eigenvalues of the transfer matrix T(u). Define A& = k 
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and A s = k s , then the free energies are defined by /&(«) 
— log k s (u). For the bulk 
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At this point we need to distinguish between the unrotated and rotated lattices. Naturally 
the bulk free energy remains unchanged in either geometry. For the rotated lattice in 
Fig. 2 the boundary face weights (A'-matrices) are gone, so the surface free energy is 
simply determined by 



k s (u)k s (u + A) 
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To obtain the surface free energy we consider the regime < u < A with < p < 1, 
where p = e~ € — > at criticality. To demonstrate the procedure, we first recall the 
derivation of the bulk free energy. We introduce the new variables 
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Up to a harmless prefactor, which we disregard in the following, the necessary conjugate 
modulus transformation of the theta function is 

0i (it, e" e ) ~ E (e" 4 ™ /e , e" 4 ^) (5.5) 

where 

CO 

E(z,x) = \{{l-x n - 1 z)(l-x n z' 1 ){l-x n ). (5.6) 

71=1 

The argument is to suppose that Kb(w) is analytic and nonzero in the annulus x < 
w < 1 and perform the Laurent expansion \ogKi,(w) = Z)^L_oo c nW n . Then inserting this 
into the logarithm of both sides of ( |5.1| ) and equating coefficients of powers of w gives 

go 



^, . ^ y-^ s i nri [2™n/e] sinh [27r(A — -u)n/e] cosh [27r(7r — 2A)n/e] (5 7) 

~y nsinh [27r 2 n/e] cosh [27rAn/e] 

Similarly, Laurent expanding log k, s (w) = Y%L-<x c n w n and solving the crossing unitarity 
relation (|5.3| ) under the same analyticity assumptions as for the bulk case gives 



^ . , ^ y-^ s i n h [47TMn/e] sinh [47r(A — u)n/e] cosh [27r(7r — 4A)n/e] (5 8) 

nsinh [27r 2 n/e] cosh [47rAn/e] 

The model with A = — ^j- in the regime < u < A and < p < 1 is the Aj, ABF RSOS 
model in regime IE II. By making use of the Poisson summation formula, the sum in 
(|5.7j) can be rewritten in terms of p. As p — > the singular part of the bulk free energy 
for L odd scales as 

f b ~p 2 - ab \ogp with a b = 2-^- (5.9) 

2\ 

with no singular contribution for L even || . 

To extract the surface specific heat exponent a s , we follow [gT], g2[ and define the local 
internal energy e s and the surface specific heat C s in the surface layer, 

df s {w,p) de s 
e s ( P ) ~ — ^— , — . (5.10) 



For the fixed boundary conditions we drop the correction energy e\ to e s (p) p^j . From 
(|5.8|) we find that as p — > 

/ s ~p 2 - as with a s = 2-^. (5.11) 

There is a logp correction factor for L odd. The surface specific heat exponent for the 
A L ABF model follows as 

a s = 2-^. (5.12) 
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Using the known scaling relation a s = a& + u, which has been explicitly confirmed for 



the eight- vertex model [IS], we expect v = (L + l)/4 for the correlation length exponent 
of the A L ABF model with odd L. 

Although the surface free energy on the original unrotated lattice depends on the 
explicit form of the boundary face weights, they play no role in the dominant critical 
singularity of the internal energy e s (p). This is seen also in the study of the eight- vertex 
model jnj. However, they do effect the correction energy e\ to e s (p). We therefore expect 
to see (presumably) height-dependent surface free energies and thus other surface critical 
exponents. 
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